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Abstract 

A number of properties of the Uehling potential are investigated. In particular, we determine 
the Fourier spatial resolution of the Uehling potential. The lowest-order correction on vacuum 
polarisation is re-written in terms of the electron density distribution function. We also discuss 
the consecutive approximations of the perturbation theory developed for the short-range Uehling 
potential in the Coulomb few-body systems (atoms). The cusp problem is formulated for few-body 
systems in which particles interact with each other by the mixed (Coulomb + Uehling) potential. 
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I. INTRODUCTION 



As is well known (see, e 



■g-, 0, Q; 



in the lowest order approximation the effect of vacuum 



polarisation between two interacting electric charges is described by the Uehling potential 
U(r) Q. In fl we have derived the closed analytical formula for the Uehling potential. For 
atomic systems this formula is written in the following three-term form (in atomic units 
h — l,m e — 1, e = 1) 
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where the notation Q stands for the electric charge of the nucleus, b = a -1 and a 
is the dimensionless fine-structure constant. Here and below h = is the reduced Planck 

Z7T 

constant (also called the Dirac constant), e is the electric charge of the positron and m e is 
the mass of the electron (= mass of the positron). 
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In Eq.([T]) K (a) is the modified Bessel function of zero order (see, e.g, 13), i.e. 
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where ip{k) is the Euler psz-function defined by Eq.(8.362) from [5]. The functions Ki\(z) 
and Ki 2 (z) in Eq.([T]) are the recursive integrals of the K (z) = K (z) function, i.e. 
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and Kir, (z) 



Ki n -x(z)dz , 



(2) 



where n > 1. 

The closed analytical formula for the Uehling potential derived in [4j allows one to obtain 
substantial acceleration of numerical algorithms used to evaluate the effects of vacuum po- 
larisation in various few-electron atoms and ions. On the other hand, the simple three-term 
formula, Eq. ([1]) , can also be used to perform theoretical analysis of the vacuum polarisation 
itself. In this study we report some interesting results of our recent investigations. These 
results can be useful in numerical applications and may lead us to better understanding of 
vacuum polarisation in atomic systems. In particular, in this study we derive the general 
expression for the lowest-order correction on the vacuum polarisation in atoms and ions writ- 
ten with the use of the electron density distribution function p(r). The arising expression 
has a number of advantages for theoretical analysis of vacuum polarisation in light atoms 



and ions. Then we determine the spatial Fourier resolution of the Uehling potential U(r). 
This leads us to the formula for the spatial Fourier resolution of the actual electrostatic field 
which is the sum of the regular Coulomb potential and short-range Uehling potential. 

Another interesting problem is the application of the perturbation theory to determine 
various corrections produced by the short-range Uehling potential. It is shown in the fourth 
Section that all consecutive approximations produced by the Uehling potential can be deter- 
mined with the use of the representations of the 0(2,l)-algebra(s). Such an algebra describes 
the dynamical representations of the original Coulomb system (or atom). We also consider 
the cusp problem for few-body systems of interacting electrically charged particles. 

II. UEHLING CORRECTION AND ELECTRON DENSITY DISTRIBUTION 

In atomic systems the lowest-order correction on vacuum polarization (AE)v.p. to the 
total energy of the LS— bound state is written in the following general form 
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(AE)v.p. = / n s (r)U(r)y L s(r)d 3 r = An / ^ LS (r)U(r)^ LS (r)r 2 dr 
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where U{r) is the central Uehling potential Eq.([T]). In Eq.([3]) the notation ^/^(r) stands 
for the bound state wave function, while ^ Ls{ r ) means the radial part of the total wave 
function. The index LS designates the angular momentum and total (electron) spin of the 
bound state. We assume that the bound state wave function ^Ls{ r ) is know. Furthermore, 
such a function can always be chosen as a real function (or it can be transformed into 
the real function). For the real wave function one finds ^Ls( r ) = ^Ls{f)- The function 
p(r) = ^* LS (r) ■ ^Ls( r ) = ^is(r) • ^/i,s(r) is the electron density distribution in the bound 
LS— atomic state. The formula which contains the electron density function p(r) can also 
be applied in those cases when the atomic system occupies one of its 'mixed' states. 

The expression for the (AE)yp, Eq.Q, can be transformed to the different form with 
the use of our formula from [4|]. Indeed, we can write Eq.([3]) as the sum of the three following 
terms Ti,T 2 and T 3 , i.e. 

(AE) V . P . = ^-(T 1+ T 2 + T 3 ) (4) 
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The explicit expressions for the Ti , T 2 and T 3 terms are reduced to the form 
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where x = 2br. The formulas for the T 2 and T 3 terms can be simplified by using a few 
different methods, e.g., integration by parts. Finally, the formula for the T 2 term takes the 
form 
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where the function is 
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The term T 3 is transformed into a similar form 
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III. FOURIER SPATIAL RESOLUTION OF THE UEHLING POTENTIAL 



In various problems from atomic and molecular physics and Quantum Elecrodynamics it is 
often important to apply the correct analytical expression for the Fourier spectral resolution 
of the Uehling potential. It is clear a priori that the Uehling potential U (r) is a static 
potential which also is a central potential. Therefore, we can expect that its representation 
as a superposition of plane waves will contain only plane waves of zero frequency. Moreover, 
it can be shown that such spatial waves are longitudinal, i.e. they are oriented along the 
spatial vector k. The Fourier spatial resolution (or Fourier resolution, for short) of the 
Uehling potential is written in the form 

d 3 k 
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where k = (k x , k y , k z ) is the wave vector and Uk(V) is the unknown spectral function. Our 
goal in this Section is to obtain the closed formula for this spectral function and compare it 
with the spectral function of the pure Coulomb potential. 

From Eq. (|12p one finds the following equation for the spectral function «k( r ) = Uk{r) 

u(k) = J J J U(r)exp(tk ■ r)d 3 r = An J j (kr)U(r)r 2 dr (13) 



where k = ^k 2 + ky + k 2 z is the radial component of the wave vector and jo(x) is the spherical 
Bessel function of zero order. Note also that the integration in the last formula is performed 
over the area occupied by the electric charges and electric field. In Eq.f JTB"]) we have used 
the fact that the Uehling potential is spherically symmetric and, therefore, we can integrate 
over all angular variables. By using Eq.flT]) and expression for the jo(x) = ^f- function one 
finds the following formula for the spectral function itfc(r) 

An r+°° 1 / 1 \ \/t 2 — 1 

where a = ^ = ^p. The integral in the last equation can be determined with the use of the 
substitution u = r^—[ - ^ reduces the integral to the form 
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Then by using the substitution y = - one transforms this integral to the form 

1 r 1 (3-y 2 )y 2 dy 5 1 2a 4 + a 2 -l, ^ n 2N 



The final formula for the spatial spectral function takes the form 
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and, since a = 4f , we find 
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As follows from the last formula the first term in the expression for u(k) is ~ k~ 3 . This gives 
the asymptotic of the spectral function at large wave numbers k. The analogous spectral 
function for the Coulomb potential between two interacting electrically charged particles is 
Uc(k) = |f (in atomic units) (see, e.g., p]). There are a number of differences between the 
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two spectral functions u(k), Eg. ([18]) . and ucir). First, the u{r) function is an approximate 
cubic function upon k^ 1 (inverse wave number), while uc{k) is an exact quadratic function. 
Second, the u(k) function, Eq. fllSp . contains inverse powers of the fine-structure constant a, 
while the uc{k) function does not depend upon a (in atomic units). These two differences 
can also be found for the corresponding electric fields which are determined as the gradients 
of the corresponding potentials. The k-component (Fourier component) of the total electric 
filed (Coulomb + Uehling potentials) is Ek = — z|f(l + U{k)) ■ k, where the function U(k) is 

and k is the wave vector. As follows from this formula each of the Ek components is oriented 
along the wave vector k, i.e. the total electric field contains only the longitudinal components 
(there are no non-zero transverse components) and these components do not depend upon 
the time t. 



IV. CONSECUTIVE APPROXIMATIONS UPON THE FINE STRUCTURE 
CONSTANT 

In general, the short-range Uehling potential can be considered as a small correction 
to the Coulomb interparticle potential. The contributions of the Uehling potential can be 
taking into account with the use of the perturbation theory. The small parameter in this 
theory is the fine structure constant a (or the factor 2 a). Based on the existence of the 
radial 0(2,l)-algebra for the original Coulomb two-body system (atom) we can find a very 
transparent physical sense for each order of approximation used in this procedure. Below 
in this Section, we restrict ourselves to the analysis of the Coulomb two-body systems. The 
cases of the Coulomb three- and four-body systems can be considered analogously. The 
Coulomb two-body system (or atom) includes the heavy nucleus with the electric charge Qe 
and electron with mass m e and electric charge — e. The corresponding Schrodinger equation 
(in atomic units) takes the following form 
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where H is the Hamiltonian, E is the energy of the system, I? is the operator of the angular 
momentum of the system (atom) and \& is the wave function. The operator L 2 depends only 



upon two angular variables (6 and 4>) and it commutes with any operator which depends 
upon the radial variable r. The explicit form of the L 2 operator is 
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The eigenfunctions of the L 2 operator are the spherical harmonics (see, e.g., Y^ m (O,0) 
and L 2 Y im (Q, <f>) =1(1 + 1)^(9, <f>). 



In 



|s| and [9 



7j (see also |8|] and [9]) it was shown that the three following operators 

1 L 2 1 L 2 

S = -r(p 2 r + — + l) , T = r Pr and U = -r(p 2 + — - l) (22) 

form the non-compact 0(2,l)-algebra with the commutation relations 

[S, T] = -%U , [T, U]=iS , [U, S] = -iT (23) 

The Casimir operator of this algebra C2 = S 2 — U 2 — T 2 equals to the L 2 operator of the 
angular momentum. The exact coincidence of the Casimir operator C2 of the radial algebra 
0(2, 1) with the Casimir operator C2 = L 2 of the algebra of three-dimensional rotations 
0(3) is related with the complementarity of the corresponding representations in the space 
of more general Sp(6N, R )- al gebra representation (so-called Moshinsky complementarity of 
representations, see, e.g., 1(3]). 

Based on this fact we can transform the operator r(H — E) = [(S + U) — E(S — U) + Q] 
from the Schrodinger equation by applying the unitary transformation exp(t/3T), where f3 
is the real parameter and T is the generator from Eq. (l22l . Indeed, by using the well known 
Hausdorff formula 

exp(-ifiT) (S ± U)exp(ifiT) = exp(±/3T) (S ± U) (24) 



Now, we can chose /3 = ln\> ' —IE and reduce the original Schrodinger equation to the form 

[v /= 2~ES-Q]0(r,e, ( />) =0 (25) 

The operator S has the discrete spectrum [ll] which is written in the general form 
S I n, £, m) = n | n, £, m), where n — 1, 2, 3, ... is an integer positive number. From Eq. (T2"5"|) 
one finds E n = — ^ (the energy spectrum in atomic units) and | n,£,m) = A<fi(r, 6,0), 
where A is an arbitrary numerical constant. All these facts and equations are very well 
known for the Coulomb two-body systems (atoms). It appears that we can modify our 



procedure for pure Coulomb systems to include the Uehling potential. Such a modification 
is discussed below. 

First, we note that from Eq. (122|) it follows that r = S + U . Second, the modified Uehling 
potential rU (r) is written in the form: 

rU(2br) = + b ^)K (2br) - jKi^br) - (^f + ~)i^ 2 (2&r)] , (26) 

This means that the Uehling potential is a function of the 2br = -(S + U) variable. Such 
a variable contain the small dimensionless parameter (the fine structure constant) a. In 
Eq.( )26|) the parameter - plays the role of the cut-off parameter, since all modified Bessel 
functions K (2br), Kii(2br) and Ki 2 (2br) decrease exponentially at large r. Briefly, this 
means that we can restrict the power series of the rU (2br) potential to a very few first terms 
and expectation values of all these terms can be determined with the use of the commutation 
relations for the 0(2,l)-algebra, Eq.( l23j) . In particular, for me it was interesting to know 
that in the lowest order approximation the effect of vacuum polarisation can be represented 
as a small change of the electric charge (screening). 

V. CUSP PROBLEM 

As is well known in any Coulomb few-body system the expectation value of the following 
operator (in atomic units) 

* = W-yi^lk (27) 

equals to the product of the corresponding electric charges q\ and q<i- In particular, for the 
electron-nuclear Coulomb interaction written in atomic units te expectation value of the 
v e N operator equals Q exactly. The operator, Eq.( l2"T|) . has a great interest for the general 
theory of bound states in Coulomb few-body systems. On the other hand, the coincidence 
of the expectation value of this operator with its expected value, e.g., with Q, is often used 
in numerical computations to test the quality of trial wave functions. 

Now, consider the case when Uehling potential is added to the Coulomb potential directly 
into the Schrodinger equation. This will lead to the change of the cusp condition between 
each pair of interacting particles. The fundamental problem is to find the = (%) value 
for few-body systems interacting by the potential which is the sum of Coulomb + Uehling 
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potential. Formally, the answer can be written in the following form (in atomic units) 



where a is the fine structure constant, 7 m 0.577215. . . is the Euler's constant. Note that 
the explicit expression for the Vij expectation value, Eq. fl28|) . is r— dependent and its finite 
limit at r — > does not exist. On the other hand, the correct definition of the cusp- value 
includes two-particle delta-function. For non-relativistic systems this delta-function is a 
constant at the distances shorter than A e = aa® = a (in atomic units), where A e is the 
Compton wavelength. In general, the non-relativistic wave function cannot produce the 
actual electron density distribution at distances shorter than A e . Instead it always gives 
the constant. Briefly, this means that we have to assume that r = aa^ = a in Eq. (l28l) . 
Therefore, the cusp expectation value (or z/y value) is finite and its numerical value is 



In other words, the Uehling potential produces a negative correction to the pure Coulomb 
cusp which is ~ a. It will be very interesting to perform highly accurate computations for 
some few-body systems with the mixed interaction potential between electrically particles 
(Coulomb + Uehling). Then we can determine the expectation values of the actual electron- 
nuclear and electron-electron cusp values. Note that in modern highly accurate computations 
the cusp values are determined to the accuracy « 1 ■ 10~ 10 — 1 ■ 10~ 12 a.u. The expected 
Uehling correction to the pure Coulomb cusp is ~ 1 ■ 10~ 3 (or 0.1 %) and it can easily be 
detected in such computations. 

VI. CONCLUSION 

We have considered a number of properties of the Uehling potential. In particular, we 
have derived the expression for the lowest-order correction to the vacuum polarisation which 
contains the electron density distribution function p(x) in atoms and/or in molecules. Such 
a representation has a number of advantages in various applications, including actual cal- 
culations of the vacuum polarisation in the mixed atomic states and variational procedures 
used to improve the electron density distribution function p(x) from 'exact' QED-based 
Hamiltonians. We also derive the explicit expression for the Fourier spatial resolution of 
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the Uehling potential. It is shown that such a resolution contains only longitudinal waves 
of zero frequency. The formula for the Fourier spatial resolution of the Uehling potential 
is investigated. Small corrections produced by the short-range Uehling potential for actual 
Coulomb systems can be taking into account by a method based on the use of the non-trivial 
0(2,l)-group for the original Coulomb few-body systems. In conclusion, we wish to note 
that the Uehling potential has many other unique properties which were never investigated 
in earlier works, e.g., the 'corrected' cusp between two interacting electric charges. Some of 
these properties will be discussed in our future studies. 
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